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Abstract 

It is shown that by making use of the Kodama vector field, as a preferred time evolution 
vector field, in spherically symmetric dynamical systems unexpected simplifications arise. In 
particular, the evolution equations relevant for the case of a massless scalar field minimally 
coupled to gravity are investigated. The simplest form of these equations in the 'canonical 
gauge' are known to possess the character of a mixed first order elliptic-hyperbolic system. The 
advantages related to the use of the Kodama vector field are two-folded although they show up 
simultaneously. First, it is found that the true degrees of freedom separate. Second, a subset 
of the field equations possessing the form of a first order symmetric hyperbolic system for these 
preferred degrees of freedom is singled out. It is also demonstrated, in the appendix, that the 
above results generalise straightforwardly to the case of a generic self-interacting scalar field. 



1 Introduction 

Let us start with a four dimensional spherically symmetric dynamical spacetime. The metric of such 
a spacetime can be given as 

ds^ ^ gABdx^dx^ + i?2(d6'2 + sine^^d^^), (1.1) 

where gAB is a non-singular Lorentzian two-metric, moreover, gAB and R are smooth functions of 
the coordinates x"^ (the capital Latin indices take the values 1 and 2) exclusively. The Kodama 
vector field (see also [3], and also for a generalisation to higher dimensional spacetimes) is 
defined as 

K" = e''''daR, (1.2) 

where e"'' denotes the volume form associated with the two- metric gAB- It can be checked by a 
straightforward calculation that 

V^i^e = 0, (1.3) 
moreover, by making use of the Einstein's equations it can also be shown (see for details |3) that 

G'fV.Kf (1.4) 
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does vanish, where Gab denotes the Einstein tensor. This later relation implies that the vector field 

ja ^ Qabj^^ 

is also divergence free. This means that even though the spacetime might be completely dynamical, 
i.e. there might be no timehke Killing vector field in the underlying setting, J° is a locally conserved 
energy fiux vector field. 



2 Dynamical spacetimes with a massless scalar field 

Let us consider now the particular case of a massless scalar field minimally coupled to gravity. The 
Einstein equations are given as 

Gab = SiTrTab (2.1) 

with energy-momentum tensor^ 

Tab - VaV^Vb^A - ^ffafcV'^VVeV'- (2-2) 

The field equation relevant for the scalar field reads as 

V*=Ve^ = 0. (2.3) 

It follows from H2.1II and 112. 2|l that 

V'^Gab = 8T:Vbi>iWe^) (2.4) 

as it should be in consequence of the diffeomorphism invariance of the underlying theory. 

Hereafter, we shall use a regular spherical coordinate system {t, r, 6*, </>) so that r is chosen to be 
the area radial coordinate, i.e. a two-sphere of radius r has proper area 

A = Anr^. (2.5) 

Then the canonical form of the line element can be given as 

ds^ = -Adt^ + Bdr^ + r^{d0^ + smO^dcj)'^), (2.6) 

where the metric functions A and B are assumed to be smooth functions of the coordinates t and 
r. Notice that in these coordinates the Kodama vector field takes the form 

K''=^(iy- (2.7) 



To recognise the intimate relation between the Einstein equations, in particular which of them 
can be derived from the others it is informative to consider the divergence 

VEab = V" {Gab - SnTab) (2.8) 



^Notice that in the particular case of the massless scalar field considered here V" {TabK'') does also vanish, as it 
should do, whenever all the independent field equations - including 12.31 - are satisfied. 
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of the tensorial expression Eab ~ Gab — STrTah vanishing of which impUes that the Einstein equations 
are satisfied. Using now the fact that Eab has only four non-zero (algebraically independent) com- 
ponents one gets that the only two non identically vanishing components - the t and r components 
- of \7°-Eab can be given as 

^2AB\ (^^^«) +2B\Eu-^^A + ABrEtr (^-^A^ +2A^Br (^-^Etr 

^ABrEtt i-g^Bj - A\Etr-g^B - A^Err-g^B + AA^B Etr = , (2.9) 

2AB^r^ i ^Efr) - B'^r^Eff^A - B^r^Efr^A + ABr^Efr ( ^b] - ABr^Err ( ^A 



-2A^Br^ ( -^Err] +2A^r^Err^B-AA^Br^Err + AA^B^Ege = Q, (2.10) 
\dr J or 

respectively. 

The standard way of interpretation goes as follows. If one assumes that the field equations 
Err = and Egg = are satisfied one can get by the substitution of these relations to Ij2.9ll and 
Il2.1fl|l a first order Hnear homogeneous system of PDEs relevant for Eu and Etr- This system 
possesses the trivial solution for trivial data which implies that the "constraints propagate", i.e. 
Ett = and Etr = are satisfied everywhere if they hold on the initial data surface. 

In the present case there is an other possible way of interpretation of the above relations, as well. 
Assume now that the field equations Eu = 0, Etr = and Err — are satisfied. The substitution 
of these relations to H2.1flll implies that the fourth Einstein equation follows from the others, i.e. 
Egg = is automatically satisfied. 

Since the other relation (I2.9|l has not been used here there has remained a further possible 
simplification between the Einstein equations Ett — 0, Etr — and Err = 0. To determine it assume 
now that only the field equations Ett = and Err = are satisfied everywhere. The substitution of 
these relations into 112.911 yields then 

ABrEtr fy^A^ + "^A^Br (y^Etr^ - A^Etr^B + AA^BEtr = 0, (2.11) 

which is a linear PDE for Etr and can be considered as ODEs for Etr on each of the time level 
surfaces. These equations possesses the trivial solution, i.e. Etr = is also satisfied everywhere, 
whenever Etr = is satisfied at the centre. 

By combining the above two statements one can conclude, as it was done by Choptuik j2| (see 
also PI), that all of the algebraically independent Einstein equations are satisfied whenever Ett = 
and Err — hold everywhere and, in addition, Etr = is guaranteed to be satisfied at the centre. 



3 The basic structure of the field equations 

By making use of the auxiliary variables 
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112. can be given in the form of the following first-order system 



5 V4(-(^n) + n(i-i?)) 

m^-- r (^-2) 

dt r ^^■^> 
Moreover, the Einstein equations Eu — 0, Etr — and E^ — can be shown to take the form 

|.B = :^il^+47rri?(n2+ $2) (3.4) 
or r ^ ' 

d , 

— B = 87rr$nvAB (3.5) 



dt 

or r ^ ' 

Because all the variables A, B, $ and 11 appears in the above equations by adapting Choptuik's 
approach one could simply evolve the scalar field according to the hyperbolic equations H3.2II and 
H3.3|l and solve the elliptic equations 113. 4|l and l|3.fi|l on each time level surface. According to the 
argument concerning the dependence of the Einstein equations on each other this method provides 
a completely satisfactory solution for the full elliptic-hyperbolic system of equations provided that 
II3.5|I is also ensured to be satisfied at the origin. This later property can easily be guaranteed by the 
obvious choice i? = 1 at the centre which also necessary to avoid the presence of a conical singularity 
there. It is a remarkable property of the gauge choice applied by Choptuik that all the equations 
113. 4|l - H3.6|l concerning the gravitational degrees of freedom are of first order ones, moreover, two 
of these equations, along with the assumption that B = 1 at the centre, determine the functions A 
and B completely. 



4 The Kodama vector field as a time evolution vector field 

The significance of the use of the Kodama vector field in relation with the above discussed dynamical 
problem gets to be transparent upon applying the Kodama vector field, as the time evolution vector 
field, instead of {d/dtY which was used, e.g., in the above described approach of Choptuik. To see 
this replace all the time derivatives with respect to the coordinate t by derivatives with respect to 
a variable r in the above given field equations according to the rule 

We can also think of the variable t as a natural parametrisation of the Kodama vector field with 
respect to which it can be given as K"^ = (d/drY . 

We also need to replace the r-derivatives in the above field equations because of the simultaneous 
change of the coordinates {t, r) into (r, p) given by the implicit relation 

T — T{t, r) and p — r. 

Then for an arbitrary function F given originally as a function of {t, r) we have the relations 

dp ' dr^ + 9i dp ' ^^■^> 
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dr dt dr ' 



where by Ij4.1|l for the coefficient dt/dr we have 

m _ 



/AB 



(4.3) 



(4.4) 



but we do not know the other coefficient dt/dp. In fact, it is not this factor what we really need. 
Instead, as it is impHed by H4.2|l and H4.3|l 



—1 

dr 



dt f dt 



d 



dp dp \dT 



d_ 
d^ 



(4.5) 



whence only the knowledge of the combined factor = ^ (^) = VAB^ is needed to get the 
field equations relevant for the the introduced new variables. This, however, also means that the 
function C will explicitly be present in the evolution equations for B, $ and 11, therefore, we shall 
need an evolution equation for C, as well. Such an equation can be derived as follows. In virtue of 
114. 4|l we have that ^/ABM^ = 1 holds. By deriving this relation with respect to p we get then that 



dp 



dt 



dr 



d 



dp 

dt_ 

d^ 



d_ 

Yp' 



dr 

AB- 



'AB 



dt_ f dt_ 
dp \dT 



\dTdp) 
1 



d_ 

d^ 



/ABj 



d_ 

d^ 



'AB 



m 

dp 



= 0, 



which in virtue of Ij4.4ll . 114.511 . (I3.4|l and (I3.6|l yields the desired evolution equation for C as 



(4.6) 



1 d 
/AB dr 



'AB = — 
2 



1 dA 
A~d^ 



1 dB 
B~d^ 



= -Anp (tf 



(4.7) 



It can also be seen that the replacements 114. l|l and 114. 5|l yields from equations 113. 2|l . H3.3|l and H3.5|l 
the relations 

^.. ^(jn) -g(i^) , BiB-i)n-c{B + i)^ ^ 



d 



dr 

d_ 
dr 



n = 



p{B2-C^) 
B{B + 1)^- C(B-l)n 



52 - C2 



p(S2 - C2) 



dr 



8 7rp<I>n. 



(4.9) 



(4.10) 



Notice that the metric variable A does not appear in either of the equations Ij4.7ll , Ij4.8ll , (j4.9ll and 
II4.1()|I . Moreover, they consist of a first order symmetric hyperbolic system for which the initial 
value problem is known to be well-posed What is even more important is that this system of 
evolution equations, along with 



^lnA-C^\nA-^-^-^-4np(^^ 
dp dr p ^ 

derived from Ij3.6ll . consist of a complete system, i.e. 



tf) =0. 



— B = 47rp (B($2 
op ^ 



n^) + 2cn $) 



B{B-1) 

p 



(4.11) 



(4.12) 



5 



derived from i^A^ and Egg — are automatically satisfied whenever 114.711 . 114. 8|l . Ij4.9|l . Ij4.1()ll and 
II4.11|I hold everywhere, moreover, II4.12|I is satisfied on the initial data surface. Notice also that 
since the equations Ij4.7ll - Ij4.1flll separates from Ij4.11ll the later one can be solved independently 
upon the functions B, C, $ and 11 are determined. 

Thus we solve the Cauchy problem relevant for the chosen dynamical system by specifying first 
initial data for the variables B, C, $ and 11 so that they satisfy the constraint Ij4.12ll . We let then the 
system to evolve under the control of the hyperbolic equations Ij4.7|l , Il4.8|l , Ij4.9ll and Il4.10|l . Finally, 
the metric variable A has to be determined by solving the Hnear PDE H4.11II in terms of the already 
known functions B, C, $ and 11. 

We would like to emphasise that the apparently singular behaviour of equations Ij4.8|l . (I4.9|l . 
Ij4.11|l and Il4.12|l at the centre, p = 0, is compensated by the fact that whenever the functions B 
and $ are smooth they can be given in a neighbourhood of the centre as 

B = l + I3p^ and $ = (4.13) 

with (3 and being smooth functions of r and p. 

It is also important to know whether the expression B^ — appearing in the denominators of 
114. 8|1 and H4.9|l may vanish anywhere. In this respect we have only the following simple observations. 
It seems to be reasonable to have initial data C = on the initial data surface which choice 
corresponds to the case where the t = const and r — const surfaces coincide at our initial data 
surface. Notice, moreover, that by 114. 7|l the value of C cannot increase during the evolution so it 
will be always non-positive. In addition, we also have by 114. 7|l and Ij4.1()|l that 

A(s + C) = -4^p(n-$)2 (4.14) 

OT 

^(B-C) =4^p(n + $)2. (4.15) 

The second of these equations implies that B — C will never vanish since S — C = 1 on the initial 
data surface. However, the first of these equation does not exclude the vanishing of B + C. It 
may happen that B + C vanishes somewhere in the computation domain. Note, however, that the 
vanishing oi B + C would correspond to a coordinate singularity of the metric given in the new 
coordinates as 

ds^ = -Idr^ - 2 — drdp + ^ ~^ dp^ + P^(d6|2 + suiO'^dct)'^). (4.16) 
B B B 

Such a break down of the new coordinate system cannot, however, occur in a neighbourhood of the 
centre where, in virtue of 114. 7|l and H4.ir)|l . B + C — 1 holds there identically. 

Finally, we would like to emphasise that to get the simple form of the evolution equation Ij4.7|l 
- I|4.1fl|l for our basic variables i?, C, $ and 11 we could also start with the metric Il4.1fi|l and apply 
the substitutions 

|v. = $ + |n, l^^'-n. (4.17) 

These later relations can be derived by making use of 113. l|l and 114. 3|l - 114. 5|l . Notice also that in this 
setup there is no need to solve any equation for A since the time evolution of all the basic variables 
is governed by H4.7II - H4.1()|l . 
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5 Final remarks 



The favourable aspects of the use of the Kodama vector field in the investigated evolutionary problem 
gets apparent if one thinks of the straightforward process yielded the above simple setup. In fact, 
it is very unlikely that whenever one is given the investigated gravity scalar field system, without 
making use of the Kodama vector field, one would ever chose the gauge in which the metric possesses 
the form Ij4.16ll . moreover, one would also apply the substitution Ij4.17ll . 

It is also important to emphasise that the above described process had immediately provided 
a system of first order symmetric hyperboHc evolution equations for the smallest possible number 
of variables. The significance of this point gets to be transparent whenever, in case of other gauge 
choices, one tries to use the standard techniques of reduction to get a system of first order PDEs from 
second order ones. It might be surprising but, in general, such a simple minded process does not au- 
tomatically yield a system of strongly hyperbolic evolution equations. Recall that without possessing 
this property even the well-posedness of the associated evolutionary problem is not guaranteed. 

It is also worth mentioning that the main conclusions of the present paper generalise straightfor- 
wardly to the case of self-interacting matter fields. In the appendix the evolution equations relevant 
for case of a generic self-interacting scalar field, along with the related modifications of the argu- 
ment appHed above for the massless scalar field can be find. The time evolution of more compHcated 
coupled gravity-matter systems, including various non-linearly coupled gauge fields like Yang-Mills 
and Higgs fields, will be studied elsewhere. 

There is a definite disadvantage of the applied canonical gauge representation, which is, in fact, 
shared by any of the gauge representations using the area radial coordinate. The disadvantage is 
that the trapped region cannot be represented in an associated coordinate domain since it necessarily 
breaks down at the marginally trapped surfaces. In particular, in case of the applied gauge choice 
the marginally trapped surfaces can only appear in the limit B = (1 — 2m /r)^^ oo. 
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Consider now the particular case of a self-interacting scalar field i.e. we shall assume that the 
dynamics of 5" can be derived from a Lagrangian of the form £ — V'^^'Ve^' + V{'^), where the 
potential V is an arbitrary but sufficiently regular expression depending only on the field variable 
For instance, V{"^) could be chosen to be the widely used quadratic potential, V{"^) = X — ^I'^) , 
with A, 5*0 > 0, although as it will be seen below the following argument does not require such a 
limitation in the functional form of V{'i>). 

The field equation and the energy-momentum tensor relevant for the above specified generic 
system read as 



Appendix 




(A.l) 



and 



Tab 



Va Wfa^- - -gab (V^ + V{^ j) . 



(A.2) 
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A direct calculation also justifies that the associated changes in equations Ij4.8|l . Ij4.9|l . Ij4.11|l and 
II4.12|I can be given as 



OT 



9 ru 9 



5 (jn) - B{B - l)n - C{B + 1)$ , AirpBjC^ ~ Bn)V + \BC{%) 

dr 



p{B^-C^) B^-C^ 
Bi-fp^) - C{fn) B{B + 1)<^-C{B-1)U 47rpB{B<f-CU)V + ^B^i^) 



B^-C^ 
d 



d 



^ Inyl-C — InA- 

Op OT 

d 

— B ^ An p (B (<i>^ 
op 



p{B^ - C2) 
(B-l) 



- An p ($- 
BV) + 2CU $) - 



-C2 

-BV)^0, 

B{B-l) 
P 



(A.3) 

(A.4) 
(A.5) 
(A.6) 



Notice first that the metric variable A does not entered to either of the formerly A-independent 
equations. In fact, the most important change is that not only the derivatives and given 
in terms of 11 and $, enter to the above field equations but also the field ^ itself via the self- 
interaction potential V^(^). This necessitates, in addition to Ij4.7|l . IIA.311 . IIA.411 and Ij4.1()ll . the use 
of the evolution equation 

for ^ . This equation can be derived from the relations defining IT and that of Ij4.1ll . moreover, \A.7\ 
is also independent of the metric variable A. More importantly, l|A.7ll . (|A.3ll . (IA.4II . (I4.in|l and H4.7|l 

can be seen to possess the form of a first order symmetric hyperbolic system 



d ^ f d 



Siu), 



(A.^ 



where the vector variable u and the symmetric coefficient matrix Ap are given as 



n 

$ 

B 



/ 





and 





c 

B 

B'^-C^ 








B 

B'-'-C'^ 
C 

B^-C^ 






\ 







0/ 



(A.9) 



moreover, the vector expression S{u) stores all the suitably arranged source terms of equations (IA.7II . 
IIA.311 . IIA.411 . H4.10|l and H4.7II . For such a system the initial value problem is known to be well-posed. 
Moreover, again we have that this system - consisting of ljA.7|l . ljA.3ll . ljA.4|l . II4.10|I and H4.7II -, along 
with ljA.5|l . comprises a complete system, i.e. ljA.6|l and Egg — can be derived from them, whenever 
ljA.7ll . ljA.3|l . ljA.4|l . II4.10|I . H4.7|l and ljA.5|l hold everywhere, moreover, ljA.6|l is satisfied on the initial 
data surface. Since, as it was before, IIA.7|l . ljA.3ll . IIA.4|l . Il4.10|l and Ij4.7|l separate from ljA.5|l the 
later can be solved for the metric variable A independently upon the functions ^f,!!, $,S and C 
have been determined. 
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